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ABSTRACT 


Intense  and  persistent  cyclonic  vortices  frequently 
dominate  the  atmospheric  circulation  to  high  levels  over 
Arctic  areas.  With  the  aid  of  satellite  infrared  images, 
this  study  examines  some  aspects  of  these  systems. 

From  analyses  prepared  by  the  Canadian  Meteorological 
Centre  (C.M.C.)  500-mb.  cyclone  tracks  are  plotted  for 
a  twelve-month  period  from  June,  1973  to  May,  197^. 

The  NOAA  ITOS  Satellite  images  used  in  this  study  were 
received  by  the  Satellite  Laboratory  at  the  University 
of  Alberta.  A  brief  discussion  of  some  characteristics 
of  these  images  is  presented,  as  well  as  some  basic 
information  on  the  recognition  of  cloud  types  and  clear 
areas. 

Selected  cases  were  chosen  from  the  summer  months 
of  July,  August  and  September,  1973.  Moisture  and 
temperature  distributions  are  derived  using  a  simple 
cloud-band  model  for  the  cloud  masses  surrounding  these 
deep  vortices. 

Two  methods  of  vertical  velocity  computation  are 
examined  in  this  study,  namely  the  vorticity  and  the 
omega  equation  techniques.  A  12  x  12  grid,  with  a  vrid- 
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length  of  381  km. ,  is  established  over  the  Canadian 
Arctic.  A  test-field  of  geopotential  gridpoint  height 
data  is  used  in  the  vertical  velocity  computations  for 
comparison  of  the  two  techniques.  C.M.C.  gridpoint 
height  data  are  then  used  to  compute  vertical  velocities 
and  compare  the  two  methods  for  a  two-week  period  in 
July,  1973*  Combinations  of  vertical  velocity  and 
moisture  ranges  are  presented  in  tabular  form  for  clear 
and  cloudy  areas  determined  from  the  appropriate  satellite 
images.  For  the  vorticity  method,  correlations  of  cloud 
with  moisture  and  ascent,  or  of  clear  areas  with  descent 
and  no  moisture  appear  poor.  From  the  omega  equation 
results,  the  ascent  -  cloud  and  descent  -  clear 
correlations  are  Improved.  Both  clear  categories  for 
moisture  thickness  (M)  0  to  2000  feet  show  more  subsidence 
than  ascent.  For  middle  cloud  with  M  £  8000  feet,  more 
cases  of  ascent  are  observed.  However,  for  low  cloud 
with  2000  £  M  ^  8000  feet,  more  subsidence  occurrs. 

Finally,  a  case  study  is  examined  to  relate  synoptic 
features  and  cloud  masses  from  satellite  images. 
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CHAPTER  I 


SYNOPTIC  AND  THEORETICAL  CONSIDERATIONS 


1.1  Introduction 

Since  500-^b.  cyclonic  systems  occur  in  all  seasons 
and  often  dominate  the  circulation  over  the  entire 
Canadian  Arctic,  any  new  information  on  the  structure 
of  these  systems  is  valuable.  The  analysis  and 
forecasting  of  meteorological  parameters  associated 
with  such  systems  is  often  complicated  and  in  some  cases 
all  but  impossible  because  of  a  lack  of  surface  and 
upper  air  data.  Fortunately,  with  the  launching  of 
meteorological  satellites  in  recent  years,  the  situation 
has  somewhat  improved,  and  satellite  photos  are  now 
available  to  augment  the  sparse  data  from  the  few 
conventional  observation  sites  in  the  Arctic. 

This  investigation  consists  of  four  main  sections. 
The  first  section  presents  monthly  charts  of  500-mb. 
cyclone  tracks  for  the  period  from  June,  1973  to  May, 
1974.  These  charts  are  chosen  to  illustrate  the 
persistence  and  frequently  erratic  motion  of  these 
systems.  The  second  section  deals  with  the  moisture 
and  temperature  distributions  in  various  areas  near 
the  upper  disturbance,  according  to  a  cloud-band 
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model  described  In  Chapter  II.  The  third  section  is 
concerned  with  the  patterns  of  vertical  velocity  and 
moisture  as  manifested  in  satellite  images  over  the 
Arctic.  The  final  section  presents  a  case  study  in 
which  the  satellite  cloud  observations  are  compared  to 
several  analyses  of  meteorological  fields  at  various 
levels. 

1.2  General  Specifications  of  Satellite  Photos  Used 

The  satellite  photographs  used  in  this  study  were 
received  and  processed  at  the  Meteorological  Satellite 
Laboratory  of  the  University  of  Alberta.  The  pictures 
obtained  were  from  the  U.S.  National  Oceanic  and 
Atmospheric  Administration  satellite  NOAA-2,  orbiting 
at  790  nautical  miles  above  the  Earth.  Since  only  infra¬ 
red  images  were  used  in  this  study,  the  following 
discussion  will  centre  on  the  characteristics  and 
specifications  of  the  Infrared  Images.  The  material 
for  this  and  the  subsequent  section  of  cloud  analysis 
has  been  taken  from  Anderson  et  al  (1973)* 

The  satellite  scanning  radiometer  measures 
radiation  in  the  10. 5  to  12.5  Jiro.  spectrum  range  with 
the  measurements  converted  to  different  shades  of  grey. 
Since  the  long-wave  radiation  emitted  by  the  Earth’s 
surface  and  clouds  varies  according  to  the  temperature 


of  the  radiating  surface,  the  darkest  areas  on  the 
infrared  photos  correspond  to  the  warmest  radiating 
surfaces  while  the  whitest  areas  correspond  to  the 
coldest  surfaces. 

The  scanning  radiometer  scans  perpendicularly 
to  the  satellite  path,  from  horizon  to  horizon.  This 
causes  considerable  distortion  of  the  picture  along  the 
edges  and  reduces  the  reliable  data  portion  of  the 
picture  to  approximately  80  percent  of  the  distance 
from  the  centre  of  the  picture  to  its  edge.  In  Figure  1 
below,  the  geometry  of  the  scan  is  shown  with  the  Earth's 
curved  surface  displayed  on  the  flat  image  plane  of 
the  infrared  scanner.  The  image  plane  is  divided  into 
eight  equal  segments,  with  the  number  above  each 
segment  indicating  the  length  in  nautical  miles  of 
the  arc  of  the  Earth's  surface  Included  in  that  segment. 
The  numbers  shown  below  the  curved  Earth's  Surface 
indicate  the  instrument  resolution.  It  ranges  from 
4.0  nautical  miles  to  14.7  nautical  miles,  beyond  which 
the  resolution  decreases  markedly  while  the  distortion 
Increases  rapidly.  This  restricts  the  reliable  data 
region  to  the  inner  860  nautical  miles  of  the  scan, 
which  is  approximately  80  percent  of  the  distance  from 
the  image  centre  to  the  image  edge,  as  stated  previously. 


Region  of  Reliable  Data 


Image  Edge 


Image  Edge 


Fig.  1.  Scan  of  the  Earth's  curved  surface  displayed 
on  the  flat  image  plane  of  the  infrared 
scanner.  (Anderson,  R.K.  et  al  1973*  p.4). 


1.3  The  Use  of  Satellite  Photos  in  Cloud  Analysis 


The  following  remarks  concern  the  appearance 
of  low  and  middle  cloud,  and  clear  areas  as  recorded 
on  infrared  satellite  images.  The  basis  for  the 
discrimination  between  low  and  middle  cloud  is  the 
different  shades  of  grey  which  cloudy  areas  show 
relative  to  cloud-free  land  or  water  areas. 

Low  cloud  areas  such  as  stratocumulus  or  stratus, 
not  obscured  by  higher  clouds,  generally  appear  grey, 
because  the  tops  of  the  low  clouds  are  cooler  than  the 
land  areas,  which  appear  as  darker  shades  of  grey. 
However,  problems  can  occur  with  this  method  of 
discriminating  low  cloud  areas  when  the  tops  of  the 
radiating  clouds  are  at  a  temperature  nearly  equal 
to  the  adjacent  land  or  water  areas,  because  the  clear 
and  cloudy  areas  exhibit  then  similar  shades  of  grey. 
This  difficulty  is  alleviated  to  some  extent  in  the 
Canadian  Arctic  by  the  large  number  of  islands  and 
areas  of  open  water  in  the  summer,  which  provide  good 
contrast  to  the  coastlines.  Thus,  an  overlying  area 
of  low  cloud  which  may  be  radiating  at  the  land  or 
water  temperature  will  still  obscure  the  coastlines, 
and  so  be  readily  differentiated. 

However,  with  winter  and  spring  satellite  images, 

the  cloud  determination  is  much  more  difficult. 
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Absence  of  a  strong  Earth-cloud  temperature  difference 
would  make  the  cloud  type  difficult  to  ascertain  from 
the  grey  shade  variations  alone.  Also,  with  only  a 
few  open  leads  present,  the  temperature  discontinuity 
along  coastlines  cannot  readily  be  used  as  a  clue  to 
the  presence  of  cloud. 

Middle  cloud  areas  such  as  altostratus,  altocumulus 
and  nimbostratus  appear  light  grey  to  white,  because  of 
their  generally  high,  cold  radiating  tops.  The  main 
problem  in  identifying  middle  cloud  is  possible  confusion 
with  cirrus  or  cirrostratus ,  since  cirriform  cloud  types 
can  also  appear  as  light  grey  or  white  images.  Erroneous 
classification  could  result  if  no  middle  cloud  layer 
existed  under  the  high  cloud.  However,  this  may  not 
be  a  serious  problem,  since  frequently  well-organized 
bands  of  thick,  high  cloud  have  layers  of  cloud  below 
them  also.  Finally,  the  determination  of  clear  areas 
on  satellite  Images  is  a  relatively  straightforward 
process,  since  the  land  areas  in  summer  are  generally 
warm,  and  consequently  appear  very  dark  in  the  photos. 

The  coastlines,  as  previously  mentioned,  also  show  a 
good  contrast,  a  fact  that  further  assists  in  determining 
the  cloudless,  areas.  Figure  2  shows  a  sample  image 
illustrating  low  and  middle  cloud,  and  clear  areas. 
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Geographical  outlines  have  been  drawn  only  where 
cloud  obscures  coastlines.  In  the  cloudless  areas, 
several  islands  are  clearly  visible,  as  for  example, 
the  northern  parts  of  Victoria  Island  (1)  and  Banks 
Island  (2).  Low-cloud  covers  southern  Victoria 
Island  (3)  and  northern  Ellesmere  Island  (4).  Here 
the  shading  of  the  cloudy  areas  is  a  light  grey, 
similar  to  that  of  the  open  water,  or  areas  between 
islands  covered  by  ice.  However,  the  absence  in  such 
areas  of  the  contrasting  coastlines  means  that  they 
are  obscured  by  cloud.  A.  middle  cloud  area  is  that 
labelled  (5)»  where  the  higher,  colder  tops  radiating 
at  a  lower  temperature  are  producing  a  shade  of  almost 
pure  white. 

1.4  Vertical  Velocity  Calculation 

1.4,1  Description  of  Flnlte-Dlf ference  Grid 

The  12  x  12  grid  used  in  the  finite-difference 
calculations  is  shown  in  Figure  3  below.  The  grid 
interval  is  381  km.  at  the  60°N  latitude  standard. 

The  inner  dashed  line  indicates  the  8x8  gridpoints 
at  which  vertical  velocities  were  calculated.  The 
upper  left  hand  corner  of  the  grid  is  labelled  (1,1) 
with  the  i  and  j  indices  increasing  along  the  directions 


of  the  arrows  shown. 
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1*1?.  3»  ^ap  of  12  x  12  arid  (outer  solid  line) 
and  8x8  grid  (inner  dashed  line). 
Oriffin  of  grid  is  labelled  “l,!” . 


1»4«2  Vertical  Velocity:  The  Vortlclty  Method: 

One  of  the  methods  used  to  compute  vertical 
velocities  is  the  vorticity  method.  The  derivation 
given  below  is  based  on  that  given  by  Petterssen  (1956). 

The  vorticity  equation  referred  to  isobaric 
surfaces  using  pressure  as  the  vertical  coordinate  is: 


U  SQp  -f  it 3Qp  4*1 U  ^Qp=  — Qp'j' -f 


1.1 


where  Op  is  the  absolute  vorticity  with  respect  to  a 
constant  pressure  surface, 

Dp  is  the  divergence  with  respect  to  a  constant 
pressure  surface, 

U)=  Jp/di  is  the  generalized  vertical  velocity 
using  pressure  as  the  vertical  coordinate, 

<^x  and  ^ y  are  the  horizontal  components  of  the 
relative  vorticity, 

and  U,  and  IT  are  the  horizontal  wind  components. 

The  equation  of  continuity  in  an  (x,y,p,t)  system  is 


Dp=-iin 

*P 

Substituting  (1.2)  into  equation  (1.1),  and  dropping 
the  subscript  p  results  in 


ip  ip  T  ix  Jr  iy 


From  the  relation 


km  -  Q2  w  2Q 

ip  V2  dp 


1.3 


-  QduJ  _  n/dQ 

dp  dp 


1.4 


It  follows  that 


are  usually  small  and  so  may  be  neglected: 

oJ.  _  uu  _  f  1 r*%y  S P  . 

Wl  Qo  o  Q* 

where  subscript  1  refers  to  an  upper  level, 
subscript  0  refers  to  a  lower  level. 

The  mean  value  theorem  may  be  used  to  evaluate 
the  Integral  as  follows: 

nil  a  Ufa -f. [  ± a +ir'd%y)  (p,-p o) 

Qi  Qo  L  QA 

where  the  bar  refers  to  the  mean  of  the  function 
within  the  square  brackets  between  the  levels  p, 
and  p0 . 

The  lowest  level  is  taken  to  be  1000  mb.  and 
the  vertical  velocity  at  1000  mb.  is  assumed  to  be 
zero.  The  values  of  u)  were  computed  at  850  mb. , 

700  mb. ,  500  mb.  and  300  mb. 

In  order  to  evaluate  equation  (1.7),  the  vorticity 
must  be  known;  the  true  vorticity  can  be  approximated 


by  calculating  the  sreostrophic  vorticity 

Q3*yV22  ++  ! 

where  V  2  is  the  Laplacian  of  the  geopotential  height 
field  H , 

f  is  the  Coriolis  parameter, 
and  g  is  the  gravitational  acceleration. 

The  wind  components  u  and  v  are  calculated  from  the 
geostrophic  wind  equations: 

The  computations  are  performed  using  finite- 
difference  formulae  as  shown  below  for  a  gridpolnt 
(l,j)  as  in  Figure  4.  The  finite  difference  formulae 
evaluated  at  gridpolnt  (i,j)  are: 


+f(l,j )  1-12 


m  t i,j ) 


1.13 
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x  i,  i+1 


Fig.  4.  Gridpoint  (i,j)  and  four  neighbouring 

gridpoints  required  for  finite-difference 
terras  in  vorticity  method  of  vertical 
velocity  calculations. 


(I^j.  =  >} ) 


1.14 


1.15 


where  HM  Is  the  geopotential  height  of  a  pressure 
level  at  point  (i,j), 

is  the  map  scale  correction 

H-Sinl  W 

factor  at  latitude  V  for  the  gridpolnt  ( 1 , j ) 
for  a  polar  stereographic  projection  at  60°N. 

is  the  Coriolis  parameter  at  gridpolnt  (i,j), 
H  is  the  grid  spacing  =  391  km. , 
is  the  time  interval  =  12  hours, 

Qj(++Ai)  is  the  geostrophlc  vorticity  at  time  "ktAt  , 
and  QsOt-ai)  is  the  geostrophlc  vorticity  at  time 

1.4.3  Vertical  Velocity:  The  Omega  Equation: 

The  second  method  to  be  used  in  this  study  for 
calculating  vertical  velocities  is  the  omega  equation 
for  a  quasi-geostrophlc  system.  Derivation  of  the 
equation  follows  that  given  by  Haltiner  (1971).  The 
quasi-geostrophlc  vorticity  equation  for  an  (x,y,p,t) 
coordinate  system  1st 


1.16 


aV'li  +V.VQ=TlMl 

v  it  »p 


where  "f  Is  the  average  Coriolis  parameter  over  the 


area  of  interest. 

The  thermodynamic  equation  assuming  no  diabatic  heat 
term  is 


a  /izy  v.  Vii  +o*ai= o 

WUpJ  ip 
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is  the  static  stability. 


To  eliminate  the  time  derivatives  in  (1.16)  and  (1.17), 
differentiate  (1,16)  partially  with  respect  to  p  and 
take  the  Laplacian  of  (1.17)  and  subtract  the  two 
resulting  equations.  The  result  of  this  manipulation 
is  the  omega  equation. 


1.18 


The  advectlon  terms  V.70  and  v.m  may  be 


written  as: 


using  the  quasi-geostrophic  wind  components 


V.VQ-q_T(2,Q) 


■ 


where  J  Is  the  Jacobian  operator. 

Similarly, 

VJiZ-J (2,&\ 

ip  \  ip) 

Thus,  the  omega  equation  may  be  written  as: 

'F‘“f 

In  order  to  solve  this  equation  numerically, 
finite-difference  forms  of  all  derivatives  in  (1.19) 
must  be  used.  The  horizontal  grid  is  identical  to 
that  used  for  the  vorticity  method,  as  shown  in 
Figure  3*  except  that,  for  the  omega  equation,  the 
dimensions  were  increased  to  20  x  20  by  adding  4 
grldpoints  to  all  sides  of  the  12  x  12  grid. 

Figure  5  shows  a  central  gridpoint  (i,j)  and  its 
eight  neighbouring  grldpoints,  as  used  in  the  horizontal 
f inite*-dif f erence  approximations.  The  two  horizontal 
operators  are  the  Laplacian  and  the  Jacobian. 

(1)  The  Laplaclan: 

The  finite-difference  form  of  the  Laplacian 
of  a  quantity  Z  evaluated  at  a  point  (1, j)  is: 

T*  2  (i,j)  *(n(i,j-i)+2(i,j+ 

where  r  refers  to  the  finite-difference  form  of  the 


Laolacian. 


,• .  X“ 


' 


X 


i-l,i 


x 


x  i-1,  j  +  1 


i,  j-1  x 


•i  I 

x 
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H 


xi,  j  +  1 


i+1 ,  j-1  x 


x 

i+U 


x  i+1,  j+1 


Pig.  5*  Gridpolnt  (i,j)  and  eight  neighbouring 

gridpoints  used  in  finite-difference  terms 
for  the  omega  equation. 


18 


(11)  The  Jacobian; 

The  general  definition  of  the  Jacobian  of  two 
variables  A  and  B  Is: 

The  finite  difference  form  of  the  Jacobian  at 
point  (i,J)  is: 


1.21 


1,4,4  Vertical  Structure  for  Omega  Equation: 

Figure  6  shows  the  isobaric  levels  required  for 
the  omega  equation  solution.  The  levels  indicated 


by  Z(l)—2fo)  are  the  levels  at  which  the  geopotential 


heights  are  used  to  solve  for  the  u/s  at  the  Intermediate 
levels  of  775  mb. ,  600  mb. ,  and  400  mb.  The  lower 
(1000  mb.)  and  the  upper  (200  mb.)  boundary  u/-  values 
are  set  to  zero. 

The  finite-difference  equation  for  the  first 
derivative  with  respect  to  p ,  for  a  quantity  iu  , 
using  the  vertical  indexing  as  in  Figure  6  evaluated 
at  say  level  2  (775  mb.)  is  the  following: 


Iim)  = 

Up4  pu)-p(3) 


1.22 


w  s  0 - 


200  mb. 


Z  (4) - 300  mb. 

w(4) - —  —  —  —  —  —  —  —  — - 400mb. 

Z(3)« - - - 500  mb. 

w(3) - — - —  —  —  — - 600mb. 

1(2) - 700  m  b. 

w(2) - 775  mb. 

Z(l) - 850  mb. 

w(l)-0 - - lOOOmb. 


Fla-.  6.  Vertical  structure  for  omega  equation  solution. 


where  the  subscript  2  refers  to  the  _&U2  evaluated 

ap 

at  level  2  (775  mb.). 

The  second  derivative  with  respect  to  p  at 
level  2  (775  mb.)  is  given  by  the  following  expression: 


1.4.5  Omega  Equation  Solution: 

The  computer  program  used  for  the  omega  equation 
solution  was  written  by  Schram  (1974)  and  modified  by 
Hume  (1975)*  Since  the  omega  equation  involves  the 
Laplacian  of  a  Jacobian  term,  the  equation  may  be 
solved  only  at  the  inner  16  x  16  gridpoints.  The 
solution  for  the  equation  was  obtained  by  a  relaxation 
procedure.  This  involves  setting  all  gridpoints  to 
a  first-guess  omega  field,  and  then  calculating  at 
each  gridpoint  a  residual  which  is  the  difference 
between  the  finite  difference  values  of  the  left  and 
right  side  of  equation  (1.19) •  If  this  residual  is 
larger  than  a  preassigned  tolerance  limit,  then  a 
correction  is  applied  to  that  gridpoint  by  adding 
to  it  a  fraction  of  the  residual.  The  entire  16  x  16 
grid  is  scanned  in  this  manner,  beginning  at  the  lowest 
level  at  which  u>  is  to  be  calculated  i.e.  775  mb.  and 


endin'?  at  the  uppermost  level,  namely  400  mb.  After 
each  complete  scan,  the  number  of  times  that  the 
residual  exceeds  the  tolerance  limit  is  recorded. 

When  this  number  is  less  than  a  prescribed  constant 
(  in  this  case  15  ) ,  the  scanning  procedure  is  halted 
and  the  current  u>  -field  becomes  the  final  solution. 

As  previously  Indicated,  the  vertical  velocities 
at  the  lower  boundary  (lOOOmb.)  and  the  upper  boundary 
(200  mb.)  were  set  to  zero.  The  lateral  boundary 
grldpoints  (the  outer  two  gridpoints  around  the  entire 
20  x  20  grid)  were  set  to  initial  first-guess  omega 
values;  these  values  were  left  unchanged  throughout 
the  subsequent  relaxation  process. 

The  static  stability  used  in  the  omega  equation 
was  calculated  from  the  following: 

<ra_LlI  M 

&  Zp  bf> 


=  R I  A& 

pe*  ap 


1.  24 


where 

equation  and  the  equation  of  state  respectively.  The 
static  stability  was  calculated  at  775  mb. ,  600  mb. , 
and  400  mb.  from  the  U.S.  Standard  atmosphere  values 
for  T  and  0-  ,  as  tabulated  by  Hess  (1959).  The 
<r  -values  at  each  pressure  level  were  left  constant 


3f 


and 


% 


from  the  hydrostatic 


. 


at  all  points  on  the  grid.  Haltiner  et  al  (1963) 
have  calculated  vertical  velocities  from  the  omega 
equation  using  three  different  methods  of  obtaining 
static  stabilities.  The  first  method  is  to  use  a 
single  static  stability,  which  is  the  average  of  all 
<r -values  calculated  at  all  grldpoints  and  at  all 
pressures  concerned.  The  second  method  involves  a 
cr -value  which  is  a  function  of  pressure  only,  and 
is  the  average  of  <r  -values  calculated  over  the  grid 
at  each  pressure  level.  The  third  method  makes  use  of 
o" -values  which  varied  horizontally  and  vertically. 
Haltiner  et  al  found  that  the  vertical  velocities 
varied  by  only  about  10  percent  in  the  lower  levels, 
but  by  larger  amounts  at  the  higher  levels.  Thus, 
the  variations  in  omega  values  due  to  the  static 
stability  variability  are  relatively  small,  so  that 
the  use  of  a  constant  static  stability,  at  each  pressure 
level,  as  calculated  from  the  U.S.  Standard  Atmosphere 
appears  justified. 

hA.  Statistical  Treatment  of  the  Results 

Several  distributions  of  moisture  and  temperature 
parameters  are  presented  in  Chapter  III.  The  statistical 
calculations  that  were  used  are  described  in  the 
following  subsections.  Brooks  et  al  (1953)  Is  the 


. 


source  for  the  following  definitions. 

(1)  Meant 


The  mean  Is  defined  as: 
n 

x  »  i/n  2  Xi 

i«  I 
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where  Xi  are  the  sample  values, 

and  11  is  the  number  of  observations  In  the  sample. 

(2)  Mode: 

The  mode  Is  defined  as  the  most  frequently 
occurring  value  of  the  variable. 

(3)  Median: 

The  median  is  defined  as  the  value  of  the  variable 
above  and  below  which  one  half  of  the  observations 
occur. 

(4)  Moments: 

The  general  definition  of  the  k+b-  moment  is 
given  by: 


yUh  »  ±  2  ( %L-X)* 

n  4*1 
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(5)  Skewness: 

The  skewness  of  a  distribution  Is  the  measure  of 
the  departure  from  symmetry  of  a  frequency  distribution. 
If  the  mode  of  the  distribution  falls  to  the  left  of  the 
mean,  the  skewness  is  positive  while  if  the  mode  falls 
to  the  right  of  the  mean,  the  skewness  is  negative. 


The  mathematical  expression  for  skewness  is: 


(6)  Kurtosis: 

Kurtosis  is  a  distortion  from  symmetry  of  the 
distribution  which  occurs  when  a  disproportionately 
large  or  small  number  of  observations  occur  in  the 
ranges  between  the  mean  and  the  extreme  values  of  the 
distribution.  In  the  former  case,  the  distribution 
is  "platykurtic"  while  in  the  latter  case  it  is 
"leptokurtic. "  The  mathematical  expression  for 
kurtosis  is: 


&  -JLl.  -  3. 

c 
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When  #i>0  ,  the  distribution  is  "leptokurtic”  while 
when  #x<0  ,  the  distribution  is  "platykurtic. " 

Figure  7  shows  schematic  distributions  illustrating 
skewness  and  kurtosis. 

(7)  "Log-Normal"  Transformation: 

Brooks  et  al  (1953)  discuss  a  method  of  transforming 
a  variable  which  is  not  normally  distributed  into  a 
variable  which  is  more  normally  distributed.  The 
transformation  is  called  "log-normal"  and  may  be 
used  when  a  distribution  has  a  positive  skew  and  is 
leptokurtic,  so  that  as  given  by  equation  (1.27) 
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Fig.  7.  Schematic  diagrams  to  Illustrate  skewness  and 
kurtosls. 
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and  as  given  by  equation  (1.2b)  are  both  positive. 

An  example  of  this  transformation  is  given  in  Appendix  A. 
The  form  of  the  transformation  is* 

U*a+bln((x-x)+C))  1.29 

where  a,  b,  o  are  constants, 

u  is  the  transformed  variable, 

x  is  the  original  variable, 

and  x  is  the  mean  of  the  original  variable. 

To  fit  the  original  distribution  of  (x-j)  ,  the  value 

of  /U-LL )  must  be  found  for  each  (x-x)  ,  where  u  is 

l  <ru  ) 

the  mean  of  u  and  is  the  standard  deviation  of  u. 

In  order  to  obtain  j  from  the  original  distribution 
of  fx-x)  ,  the  results  of  Finney  (19^1)  may  be  used. 

I U-U  \=  1  f  lh / I-ku  (x-xj)  -ht/  7 

I  (Tu  /  L  ’  /  £  ml 

iu=p.  y,  +  JjX+iJl/3  +f±  X,  -|rr+i 

If)  (l+U>*)  1.32 

where  6%  is  the  standard  deviation  of  the  original 
distribution  of  x, 

ft  is  the  skew  of  the  original  distribution 
of  x. 

The  procedure  Involved  to  obtain  the  transformation 
of  x  to  is  the  following. 


r 


. 


(a)  Calculate  &  and  Zf*  to  determine  whether  they 


are  both  positive  and  hence  If  the  " log-normal" 
transformation  may  be  applied. 

(b)  Calculate  an  and  for  the  distribution  of  x. 

(c)  Calculate  <ru  from  equation  (1.32). 

(d)  Calculate  ^u-u  j  from  equation  (1.30). 

Using  the  values  of  j  obtained  by  the  above 
procedure,  a  table  of  probability  integrals  for  the 
normal  frequency  curve  such  as  Appendix  II  (c)  of 
Brooks  et  al  (1953)  ®ay  be  used  to  obtain  the  frequency 
distribution  of  W)  • 

"Log-normal”  distributions  are  particularly  useful 
in  statistical  treatment  of  daily  rainfall  data. 


CHAPTER  II 


SOURCES  AMD  METHODS  QF  OBSERVATIONS 

The  sources  of  data  and  the  manipulations  required 
to  process  the  data  for  the  first  three  parts  of  this 
study  are  considered  below.  The  maps  required  for  the 
case  study  are  described  in  Chapter  IV. 

2,1  500-mb.  Cyclone  Tracks 

The  first  objective  of  this  study  was  the  determination 
of  the  monthly  500-mb.  cyclone  tracks  over  the  Arctic 
areas  for  a  period  of  one  year,  from  June,  1973  to 
May,  1974.  The  charts  used  for  this  purpose  were  the 
500-mb.  height  and  thickness  analyses  issued  by  the 
Canadian  Meteorological  Centre  (C.M.C.)  dally  at  0000  Z 
and  1200  Z  (Greenwich  Mean  Time).  The  criteria  used 
when  recording  the  track  of  a  500-mb.  low  are  the 
following. 

(I)  The  low  centre  was  labelled  as  such  on  the  C.M.C. 

500-mb.  charts. 

(II)  The  low  persisted  for  at  least  three  consecutive 
twelve-hour  periods. 

(III)  The  centre  of  the  low  was  north  of  60°N,  and 
between  approximately  40°W  and  180°W,  longitude.  Over 
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the  Arctic  Ocean,  this  rule  was  modified  to  the  extent 
that  the  centre,  in  order  to  qualify,  needed  to  be 
within  350  to  400  nautical  miles  of  the  nearest 
radiosonde  station. 

The  500-mb.  cyclone  tracks  for  this  period  are 
shown  in  Figures  8  to  13.  Since  the  period  of  study 
is  only  a  single  year,  generalizations  concerning  the 
seasonal  or  monthly  variations  in  position,  frequency 
and  persistence  of  these  systems  are  not  warranted. 
However,  a  few  specific  observations  concerning  the 
year  under  study  may  be  made.  Although  the  overall 
motion  of  the  majority  of  cyclones  was  from  west  to 
east,  the  paths  followed  by  individual  systems  frequently 
appeared  to  be  erratic  or  at  least  complex.  This 
complexity  may  in  part  sitnply  reflect  the  sparseness 
of  data.  Thus,  the  path  of  a  low  may  appear  to  loop, 
while  the  centre  has  actually  moved  in  a  relatively 
uniform  manner.  An  example  of  such  a  questionable 
looping  movement  occurred  between  02/1200  Z  and 
03/1200  Z,  June,  1973  in  the  vicinity  of  78*N,  90°W. 

Since  this  happened  over  such  a  short  period  and  the 
loop  was  very  small,  it  is  difficult  to  say  whether 
looping  is  real  or  fictitious,  l.e.  a  consequence  of 
the  Inherent  uncertainties  in  the  500-mb.  analysis. 
However,  in  some  cases,  looping  or  erratic  motion  is 


so  persistent  that  the  complicated  trajectory  is 
probably  genuine.  An  example  is  the  low  of  13/1200  Z 
to  31/1200  Z,  April,  197^  whose  path  underwent  several 
reversals  in  direction  within  a  period  of  over  two  weeks. 

It  came  as  a  surprise  to  find  that  so  few  tracks 
affected  the  Yukon-Mackenzie  River-  Great  Slave  Lake- 
Great  Bear  Lake  areas  in  the  12-month  period  analysed. 

To  determine  whether  this  is  typical  would  require 
examination  of  at  least  ten  more  years  of  data. 

2.2  Selected  Cases  of  500-mb.  Cyclone  Tracks  for 

Moisture  and  Temperature  Analysis 

For  the  second  part  of  the  study,  500-mb.  cyclone 
tracks  were  chosen  from  the  months  of  July,  August  and 
September,  1973*  The  lows  chosen  were  required  to 
have  500-mb.  positions  relatively  close  to  upper-air 
observing  stations  in  order  that  observations  of  cloud 
associated  with  the  disturbance  would  be  available  as 
well.  Another  requirement  was  that  the  500-mb.  system 
selected  for  study  was  the  dominant  cyclone  disturbance 
affecting  the  Canadian  Arctic  for  the  duration  of  a 
specific  case  study.  Table  I  below  lists  the  dates  of 
500-mb.  cyclone  tracks  chosen  for  further  study. 

Plots  of  these  tracks  were  previously  presented  in 
Figures  8  and  9» 


Table  I 


500  mb.  low  cases  chosen  for  moisture  and 
temperature  analysis  using  satellite  photographs 


Case  Number 


Date,  Time  G.M.T. 


1 

2 

3 

4 

5 

6 

7 

8 


06/0000Z-07/0000Z,  July 
09/1200Z-14/0000Z,  July 
19/OOOOZ,  July-02/1200Z,  August 
10/0000Z-11/0000Z,  August 
13/1200Z,  August 
18/0000Z-24/1200Z,  August 
13/1200Z-17/1200Z,  September 
26/1200Z-29/0000Z,  September 


* 


Pig.  8.  500-mb.  cyclone  tracks  for  June,  July  1973. 
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Several  weak  lows  occurred.  The  strongest  The  dominant  low  was  17/1200  Z  to  31/1200 
was  13/0000Z  to  27/1200Z.  to  be  studied  In  detail  In  Chapter  III  an 

In  the  case  study  of  Chapter  IV. 


Fig.  9.  500-mb.  cyclone  tracks  for  August,  September,  1973. 
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The  strongest  low  was  15/0000  Z  to  24/1200  Z.  Few,  weak  systems  occurred. 


Fig.  10.  500-mb.  cyclone  tracks  for  October,  November,  1973. 
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Fiff.  11.  500-rob«  cyclone  tracks  for  December,  1973. 

January,  197*+. 
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Fig.  12.  500-mb.  cyclone  tracks  for  February,  March,  19?4. 
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Fig,  13.  500-mb.  cyclone  tracks  for  April,  May,  1974. 
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Few,  but  intense,  500-mb,  cyclones  occurred  with  the 
strongest  and  most  erratic  being  from  13/1200  Z, 
April ]  1974  to  07/1200  Z,  May,  1974. 


The  cloud  band  model,  used  In  the  analysis  of 
the  temperature  and  moisture  distributions  of  the  cases 
listed  in  Table  X,  was  adapted  from  a  model  devised 
by  Shenk  and  Brooks  (1965).  They  investigated  the 
distribution  of  several  cloud  parameters  associated 
with  extratropical  vortices  in  all  seasons,  and  for 
three  classifications  of  vortex  development:  developing, 
mature  and  dissipating.  In  the  present  study,  no 
seasonal  or  vortex  development  separations  are  made. 

The  model  used  is  shown  in  Figure  14  below.  The 
various  regions  shown  In  the  diagram  are  defined  and 
labelled  with  letters,  as  follows: 

(a)  Area  Ci 

This  classification  includes  all  cloud  areas 
around  the  vortex  centre.  The  areas  may  or  may  not 
be  organized  in  bands.  For  each  cloud  observation, 
a  further  classification  as  to  low,  C(L)  or  middle, 

C(M )  cloud  types  was  made  according  to  the  criteria 
set  out  in  Chapter  I.  The  limitation  on  the  cloud 
associated  with  a  given  vortex  centre  was  subjectively 
set  as  the  edge  of  the  cyclonic  circulation  around 
the  500-mb.  low  centre. 

(b)  Area  An: 

This  is  the  clear  area  located  to  the  north  of 
■the  vortex  centre  and  north  of  the  cloudmass  C. 


. 


J 


An 


Fig.  14.  Idealized  cloud  band  model  showing  regions 
around  vortex  centre.  Adapted  from  Shenk 
and  Brooks  (1965). 


(c)  Area  As: 


This  is  the  clear  area  to  the  south  of  the  vortex 
centre  and  also  to  the  east  of  the  cloudmass  C. 

(d)  Area  B: 

This  comprises  the  clear  area,  to  the  west  of 
the  cloudmass  or  between  the  cloud  bands.  It  must 
lie  to  the  south  of  the  vortex  centre  if  located  to 
the  west  of  the  cloud  mass.  The  exception  to  this  is 
that  if  a  cloudless  area  extends  into  the  centre  of 
the  vortex  and  continues  north  of  the  centre,  this 
area  is  then  also  classified  a  "3" .  For  the  clear 
area  between  cloud  bands,  the  "B"  classification 
applies  in  all  cases,  whether  north  or  south  of  the 
vortex  centre. 

Figure  15  below  shows  a  satellite  image  with 
several  cloud  areas  classified  according  to  the 
Shenk  and  Brooks  model.  A  main  cloud  vortex  centre 
appears  near  (1)  in  Eaffln  Bay,  with  a  secondary 
centre  near  (2)  over  central  Baffin  Island.  A  double 
low-pressure  centre  is  also  evident  on  the  surface 
analysis  (not  shown  here)  with  cloud  vortex  (1)  being 
close  to  the  main  surface  low  position  and  vortex  (2) 
close  to  a  weaker  surface  trough  over  Baffin  Island. 
The  500-mb.  chart,  analyzed  at  25/1200  Z,  July  (about 
four  hours  before  the  satellite  image  time)  shows  a 
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500-mb.  vortex  centre  at  (3)  In  Figure  15  with  a 
large  cyclonic  circulation  covering  the  entire  central 
and  eastern  Arctic  regions.  Following  the  cloud  model, 
several  areas  are  labelled  on  the  satellite  image. 

These  arei  (4)  over  western  Greenland  -  middle  cloud, 

C(M),  around  vortex  centre*  (5)  over  northern 
Greenland-clear  area.  An,  to  the  north  of  the  vortex 
centre  and  north  of  the  cloudmass  area  (4)*  (6)  over 

southern  Baffin  Island  -  low  cloud  area,  C(L),  around 

vortex*  (?)  over  northern  Hudson  Bay  -  clear  area,  B, 
to  the  south  of  the  vortex  centre,  and  west  of  the 
cloudmass  over  northern  Quebec. 

In  order  to  position  each  radiosonde  station 
according  to  this  model,  an  adjustment  for  the  time 
between  the  cloud  observation  time  and  the  radiosonde 
observation  time  had  to  be  made.  The  reason  for  this 
is  the  following:  at  the  time  of  the  satellite  pass, 
a  cloud  mass  could  have  been  present  over  the  radiosonde 
station,  so  that  the  satellite  image  would  show  the 
station  obscured  by  cloud.  However,  at  the  time  of 
the  upper-air  observation,  the  station  could  have 
been  in  a  clear  area  so  that  the  conventional  meteorological 
parameters  recorded  at  the  station  would  differ  from  those 
inferred  from  the  satellite  image.  Corrections  were  made 
to  compensate  for  this  motion  of  the  cloud  area  between 
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the  time  of  the  satellite  photo  and  the  upper-air 
observation  time.  The  scheme  used  was  to  compute 
the  speed  and  direction  of  motion  of  the  500-mb.  low 
centre  at  the  time  when  the  satellite  image  was  taken. 

This  was  done  by  using  the  low-centre  positions  at  the 
nearest  synoptic  time  periods  before  and  after  the 
photo  time.  Figure  16  shows  a  sketch  of  the  adjustment 
procedure.  In  this  example,  the  satellite  image  was 
taken  at  12/2100  Z  and  the  radiosonde  observation  was 
taken  at  13/0000  Z.  The  speed  and  direction  of  motion 
of  the  500-mb.  low  was  determined  by  using  the  low-centre 
positions  at  synoptic  times  nearest  the  satellite  photo 
observation  time  (here  12/1200  Z  and  13/  0000  Z).  Then, 
the  cloud-edge  position  was  corrected  by  moving  it  in 
the  direction  of  motion  of  the  vortex  centre  by  a 
distance  equal  to  that  covered  by  the  vortex  in  the 
three-hour  period  between  the  satellite  photo  time  and 
the  time  of  the  upper  air  observation.  The  new,  corrected 
cloud-edge  position  is  then  the  cloud  position  at  the 
radiosonde  observation  time.  Since  the  radiosonde 
observations  are  taken  at  twelve-hour  intervals,  this 
correction  did  not  have  to  be  applied  for  periods  in 
excess  af  six  hours. 

On  the  basis  of  the  500-mb.  low  tracks  recorded  in 
Table  I,  43  infrared  Images  were  chosen  for  analysis. 

A  total  of  216  radiosonde  station  soundings  were  classified 


Fig.  16.  Adjustment  of  cloud  pattern  for  satellite 

photo-upper  air  observation  time  difference. 
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into  low  or  middle  cloud  or  clear  areas,  according  to 
the  cloud-band  model  described. 

Dewpoint  depressions  were  obtained  from  the  C.M.C. 
surface,  850-mb.  700-mb.  and  500-nib.  analyses.  Vertical 
temperature  differences  between  the  surface  and  850  mb. , 

#50  and  700  mb.  and  700  and  500  mb.  were  also  recorded. 

A  third  field,  the  moisture  thickness,  was  obtained 
from  radiosonde  data  plotted  routinely  at  the  Arctic 
Weather  Central,  Edmonton  Alberta.  The  moisture  thickness 
at  a  given  radiosonde  station  was  calculated  by  totalling 
all  layers,  below  approximately  450  mb.,  whose  dewpoint 
depression  was  less  than  5aC.  The  actual  layer  thicknesses 
were  recorded  using  the  tephigram  pressure-level  heights 
directly.  It  is  realized  that  use  of  the  tephigram  to 
obtain  the  thickness  in  this  manner  introduces  an  error, 
since  the  pressure-level  heights  recorded  on  the  tephigram 
are  based  on  the  U.S.  Standard  Atmosphere.  An  estimate 
of  the  error  introduced  by  this  was  made  by  considering 
a  1000  to  500-mb.  thickness  of  5460  geopotential  metres 
for  the  U.S.  Standard  Atmosphere.  Examination  of  several 
of  the  C.M.C.  thickness  charts  Indicated  a  minimum 
thickness  of  about  5040  gpm.  and  a  maximum  thickness 
of  about  5760  gpm.  Thus,  using  the  U.S.  Standard 
Atmosphere  thicknesses,  moisture-layer  thicknesses 
near  the  deepest  500-mb.  lows  are  a  maximum  of  about 


8  percent  too  high,  and  a  maximum  of  about  5  percent 
too  low  near  the  strongest  500-mb.  ridges.  These  are 
maximum  errors;  the  actual  percentage  errors  are 
generally  less. 

Another  error  arises  in  estimating  the  thicknesses 
from  the  tephigrams  when  determining  the  limits  of 
layers  which  have  dewpoint  depressions  less  than  5°G. 
Linear  interpolation  between  the  plotted  dewpoint 
depressions  is  used  to  obtain  the  intermediate  level  at 
which  the  dewpoint  depression  is  5°C.  The  error  in  the 
layer  thickness  estimation  caused  by  this  interpolation 
is  approximately  -  lOmb.  or  -  300  feet,  depending  on  the 
pressure  level  at  which  the  interpolation  is  done. 

A  fourth  field,  extracted  from  radiosonde  data, 
was  the  minimum  dewpoint  depression.  It  is  the  lowest 
dewpoint  depression  of  the  tephigram  plot  up  to  about 
20,000  feet  above  m.s.l.,  excluding  the  surface  dewpoint 
depression. 

2.3  Vertical  Velocity  Calculations 

From  the  several  cases  considered  in  the  previous 
analysis  of  moisture  and  temperature,  the  case  of 
19/0000  Z  to  31/1200  Z,  July,  1973  was  chosen  for  a 
more  detailed  vertical  velocity,  moisture  and  cloud 
pattern  analysis. 


For  the  vertical  velocity  calculations  using  the 
vorticity  method,  the  geopotential  heights  at  1000, 

850,  700  and  300  mb,  were  required.  The  omega  equation 
vertical  velocities  required  only  the  height  data  at 
850,  700,  500  and  300  mb.  These  data  were  obtained 
from  a  magnetic  tape  of  gridpolnt  heights  supplied  by 
the  Canadian  Meteorological  Centre,  Montreal,  Quebec. 
Since  the  grid  used  in  this  study  is  identical  in 
orientation  to  that  used  by  C.M.C.,  the  gridpolnt 
data  from  the  tape  were  used  directly  as  input  for 
the  vertical  velocity  programs. 
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CHAPTER  III 


MOISTURE,  TEMPERATURE  AND  VERTICAL  VELOCITY  RESULTS 


3*1  Moisture  and  Temperature  Distributions  Associated 

With  500-^~b.  Lows  In  Simmer 

3.1*1  Dewpoint  Depressions 

The  dewpoint  depression  distributions  are  shown 

in  the  histograms  of  Figures  1?  to  20.  The  mean,  mode, 

median  and  the  number  of  observations  for  each  distribution 

are  shown  in  Table  II.  The  abscissa  labels  refer  to  the 

0 

dewpoint  depressions  such  that  e.g.  a  bar  at  " 9  " 
represents  the  percentage  of  observations  with  a 
dewpoint  depression  of  8°  or  9°C.  The  distributions 
are  also  labelled  according  to  the  cloud  model  described 
in  Section  2.2.  The  C(L)  and  C(M)  refer  to  low  and 
middle  cloud  areas,  respectively,  while  the  An,  As,  B 
distribution  denotes  the  combined  clear  areas. 

3.1.2  Surface  Dewpoint  Depressions 

Both  the  C(L)  and  C(M)  areas  have  strong,  positively- 
skew  distributions,  with  the  former  having  a  slightly 
smaller  mean.  This  appears  reasonable, since  some  of 
the  cloud  areas  could  extend  right  to  the  surface. 

Moreover,  precipitation  falling  from  the  low  cloud 
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would  also  cause  smaller  surface  dewpoint  depressions. 
The  Individual  clear-area  distributions  should  be 
considered  as  tentative,  because  of  the  small  sample 
sizes.  The  combined  (An,  As,  E)  clear  area  -  a  larger 
sample  -  has  a  considerable  drier  distribution  (mean 
dewpoint  depression:  6*C. )  than  the  cloud  categories. 
Since  1200  Z  observations  are  included,  some  of  the 
small  surface  dewpoint  depressions  in  the  clear  areas 
may  be  caused  by  nocturnal  radiation  cooling. 

3.1.3  850-mb.  Dewpoint  Depressions 

Again,  both  the  C(L)  and  C(M)  distributions  are 
more  moist  than  the  clear-area  distributions.  However, 
the  C(L)  distribution  is  slightly  drier  than  the  C(M) 
distribution.  This  is  contrary  to  what  one  would  expect, 
since  the  850-mb.  level  would  normally  be  considered 
most  representative  of  low-cloud  moisture.  It  is  likely 
that  some  of  the  low  cloud  observations  are  fog  or  low 
stratus  associated  with  the  numerous  open  water  areas 
of  the  summer  Arctic.  On  such  occasions,  the  cloud 
may  not  extend  to  the  850-mb.  level  and  the  dewpoint 
depressions  at  850-mb.  will  be  larger.  The  middle 
cloud  moist  distribution  may  be  caused  by  cloud  layers 
extending  downward  to  the  850-mb.  level.  The  clear  areas 
have  dry  distributions  with  a  mean  dewpoint  depression  of 
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Tables  II 

Dewpoint  depression  statistics 


for  low, 

middle 

cloud  and 

clear  areas 

Class 

Mean 

Mode 

Median 

No.  of  observations 

Surface 

C(L) 

2.9 

1.0 

1.8 

98 

C(  M) 

3.4 

1.0 

1.8 

72 

An ,  As , 

B 

6.0 

4.0 

4.0 

46 

850  mb. 

C(L) 

6.3 

0.0 

5.4 

97 

C(M) 

4.5 

0.0 

1.9 

72 

An,  As, 

B 

11.8 

12.0 

11.3 

46 

700  mb. 


C(L) 

9.5 

7.0 

8.0 

98 

C(M) 

5.0 

1.0 

2.2 

71 

An ,  As ,  B 

13.1 

10.0 

11.4 

45 

500  mb. 

C(L) 

11.0 

13.0 

10.7 

98 

C(M ) 

8.6 

5.0 

7.1 

71 

An ,  As ,  B 

14.0 

13.0 

12.8 

46 

Percentage  of  Cases 
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17*  Surface  dewpoint  depression  distributions. 
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Fig.  18.  850-mb.  dewpoint  depression  distributions 
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Fi*.  19.  700-mb.  dewpoint  depression  distributions 
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Fi*.  20.  500-mb.  dewpoint  depression  distributions 


3.1.4  700-mb.  Dewpoint  Depressions 

The  700-mb.  C(L)  histogram  shows  widely  distributed 
moisture.  The  observed  pattern  may  be  caused  by  moisture 
extending  above  the  low  cloud  to  the  700-mb.  level. 

The  saturated  cases  observed  at  700-mb.  in  the  C(L) 
areas  could  be  simply  the  result  of  incorrect  cloud 
type  classification.  Occasionally,  a  satellite  image 
grey-shade  was  intermediate  between  the  bright  white 
of  middle  cloud  and  the  darker  shade  of  low  cloud. 

Such  a  case  could  be  wrongly  classified  as  low  instead 
of  middle  cloud.  The  moist  C(M)  distribution  is  to  be 
expected,  since  700-mb.  is  a  representative  level  for 
middle-cloud  analysis.  Some  72  percent  of  the  cases  in 
the  C(M)  areas  have  dewpoint  depressions  6  5°C. 

The  clear  distributions  are  again  considerably  drier, 

o 

with  a  combined  An,  As  and  B  mean  of  13*1  0. 

3.1.5  500-mb.  Dewpoint  Depressions 

The  500-mb.  level  is  frequently  above  the  cloud 
layers.  The  C(M)  area  shows  a  significant  percentage 
of  cases  with  depressions  ^  5°C. ,  probably  caused  by 
middle  cloud  extending  to  this  level.  The  C(L)  and 
clear  areas  are  considerably  drier. 


3.1.6  Minimum  Dewpoint  Depressions. 


Figure  21  shows  the  minimum  dewpoint  depression 
distributions  with  the  mean,  mode  and  median  in  Table  III. 
These  results,  obtained  from  tephigrams,  are  the  minimum 
dewpoint  depressions  to  approximately  20,000  feet  M.S.L. 
The  surface  dewpoint  depression  was  excluded  since  it 
may  be  unrepresentatively  low  for  clear  areas,  at  the 
1200  Z  observation,  because  of  radiational  cooling. 

The  abscissa  labels  again  Indicate  ranges  of  minimum 

o 

dewpoint  depressions,  e.g.  a  bar  at  2.5  C.  denotes  the 
percentage  of  cases  with  minimum  dewpoint  depression 
between  2.0°C.  and 

The  histograms  show  most  cases  in  the  C(L)  and 
C(M)  distributions  in  the  0.5°C.  range.  A  considerable 
spread  of  minimum  dewpoint  depressions  occurs  with  the 
clear  areas. 

3.1.7  Moist  Layer  Thicknesses 

Moist-layer  thickness  distributions  are  shown  in 
Figure  22 j  the  corresponding  means,  modes  and  medians 
are  given  in  Table  III.  Along  the  abscissa,  a  bar  at 
»» y indicates  the  percentage  of  cases  with  a  moist-layer 
thickness  from  6000  to  8000  feet. 

The  majority  of  the  C(L)  cases  are  in  the  3#5»7 
ranges  but  some  large  moisture  thicknesses  also  occur. 
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Table  III 


Minimum  dewpoint  depression  and  moist  layer  statistics 
for  low,  middle  cloud  and  clear  areas 


Class 

Mean 

Mode 

Median  No 

.  of  observtlons 

Minimum 

dewpoint  depression  (  C.  ) 

C(L) 

0.8 

0.0-0. 5 

~0.4 

95 

C(M) 

0.6 

0.0-0. 5 

—'0.4 

64 

An ,  As ,  B 

4.8 

4. 5-5.0 

-4.7 

42 

Moist  layer  thickness 

( feet ) 

C(L) 

6695 

o  o 
o  o 
o  o 

i 

/V  4950 

95 

C(M ) 

11846 

14000 

-15000 

19000 

-20000 

- 13500 

64 

An ,  As ,  B 

1332 

0-1000 

-  750 

42 
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Pl^.  22.  Moist-layer  thickness  distributions 


Large-thickness  cases  are  difficult  to  explain,  since 
this  would  mean  that  some  grey  areas  on  the  satellite 
image  purporting  to  indicate  low  cloud  have,  on  the 
basis  of  corresponding  radiosonde  data,  much  thicker 
moisture  layers  than  expected.  This  could  be  due  to 
errors  in  identifying  the  cloud  types  from  the  satellite 
images.  However,  since  22  percent  of  the  cases  have 
moist  layers  exceeding  10,000  feet  in  thickness,  it  may 
be  partly  the  result  of  cloud-thickness  over-estimation 
when  using  a  dewpoint  depression  ^  5°C.  as  the  criterion 
for  cloud  layers.  The  majority  of  the  C(M)  distributions 
has  large  moisture  thicknesses.  The  few  small  thicknesses 
could  be  associated  with  thin  layers  of  middle  cloud. 

The  vast  majority  of  the  clear  observations  occur  in 
range  1,  as  may  be  expected. 

3.1.8  Vertical  Temperature  Differences 

The  distributions  of  vertical  temperature  differences 
are  shown  in  Figures  23  to  25.  and  their  means,  modes 
and  medians  are  listed  in  Table  IV.  The  labelling  is 
the  same  as  that  used  for  the  dewpoint  depressions. 

Few  differences  in  the  distributions  are  apparent  for  a 
given  pair  of  levels.  This  means  that  little  significant 
difference  is  observed  in  the  stability  between  pairs 
of  levels  for  the  various  cloud  and  clear  categories. 


Vertical  Temperature  Difference (°C. ) 

Fig.  23.  Temperature  (surface) 
distributions. 
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Vertical  Temperature  Difference  (°C.) 
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Fig*.  24.  Temperature  (850  mb.)  -  Temperature  (700  mb.) 
distributions. 
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Vertical  Temperature  Difference  (°C.) 

Pig.  25.  Temperature  (700  mb.)  -  Temperature  (500  mb.) 
distributions. 


Table  IV 


Vertical  temperature  difference  statistics 


for  low, 

middle 

cloud  and  clear  areas 

Glass 

Mean 

Mode 

Median  No.  < 

of  observations 

Temperature  ( Surface ) -Temperature 

oc 

o 

mb.  ) 

C(L) 

4.8 

4.5 

—  4.6 

98 

C(M) 

5.6 

5.0 

—  4.9 

73 

An ,  As , 

B 

4.8 

3.0 

-3.7 

46 

Temperature  (850  mb. ) -Temperature 

(?00 

mb.  ) 

C(L) 

6.7 

7.0 

-  6.3 

98 

C(M) 

7.5 

8.0 

-7.0 

71 

An,  As, 

B 

8.4 

7.0 

-7.7 

46 

Temperature  (700  mb. ) -Temperature 

o 

o 

mb.  ) 

C(L) 

15.4 

14.0 

-14.1 

98 

C(  M ) 

15.4 

16.0 

-15.0 

72 

An,  As, 

B 

15.8 

16.0 

-16.4 

46 

3.1.9  The  "Log-Formal"  Transformation 


The  " log-normal"  transformation  may  be  applied  to 
a  distribution  which  is  positively  skew  (V,  >o)  and 
leptokurtic  (i*>o)  in  an  effort  to  obtain  a  new,  more 
closely  normal,  distribution. 

One  transformation  will  be  done  to  illustrate  the 
procedure.  The  example  to  be  transformed  is  the  surface 
dewpoint  depression  distribution  for  the  C(L)  area  shown 
in  Figure  1?.  The  details  of  the  calculations  for  the 
transformation  are  given  in  Appendix  A.  The  resulting 
transformed  distribution  is  shown  in  Figure  26.  This 
new  distribution  appears  closer  to  normal  and  could  be 
compared  statistically  with  other,  similarly  transformed, 
distributions. 

3.2  The  Vertical  Velocity  Calculations 

3*2.1  Vertical-Velocity  Program  Test 

To  compare  vertical  velocities  calculated  using 
the  vorticity  and  omega  equation  methods,  both  programs 
were  run  on  an  identical  12  x  12  grid  of  geopotential 
heights  as  a  test  field.  Referring  to  Figure  3»  with 
the  top  edge  of  the  grid  labelled  1*1,  the  heights 
increase  uniformly  with  increasing  index  I.  On  the 
schematic  cross-section,  Figure  27#  the  geopotential 
height  gradient  increases  with  height  from  20  geopotential 
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*1*.  26.  Dewpoint  depression  frequencies  (upper) 

and  transformed  variable  frequencies  (lower) 


metres  (gpm)  per  grid  length  at  850  mb.  to  50  gpm  per 
grid  length  at  300  mb.  Thus,  the  geostrophic  wind 
speed  and  direction  at  a  given  level  are  the  same, 
but  the  speed  Increases  with  height.  The  grldpolnt 
vorticlty  differences  are  due  only  to  variations  in 
latitude.  Under  these  conditions  the  vertical  velocities 
result  solely  from  the  advection  by  the  geostrophic 
wind  of  the  Coriolis  component  of  absolute  vorticlty. 

In  Figures  28  to  31»  the  results  of  the  computations 
using  this  test  field  are  shown.  For  the  vorticlty 
method,  the  calculated  vertical  velocities  increase 
steadily  from  low  values  at  850  mb.  to  high  values  at 
300  mb, j  a  consequence  of  the  increasing  geopotential 
height  gradient  with  height.  However,  the  omega  equation 
results  show  a  slight  Increase  in  vertical  velocities 
between  775  mb.  and  600  mb. ,. followed  by  a  decrease  at 
400  mb.  A  possible  cause  of  this  discrepancy  is  that  the 
vertical  velocity  at  the  200-mb.  upper  boundary  is  set 
to  zero  for  the  omega  equation  calculations.  This 
results  in  a  large  second  derivative  of  or  with  respect 
to  p,  a  situation  which  is  remedied  in  the  solution 
of  the  omega  equation  by  reducing  the  u> -value  at  400  mb. 

At  the  lowest  computation  levels,  the  vertical 
velocities  obtained  from  the  two  methods  are,  however, 
in  excellent  agreement.  For  example,  at  I  =  6,  J-  6, 
the  775-rab.  omega  equation  velocity  is  0.031  jib.  sec., 
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the  same  value  as  the  average  of  the  700-  and  850-mb. 
vortlcity  method  vertical  velocities.  For  this  study, 
the  7?5-  and  600-mb.  u>  -  values  are  used  so  that  errors 
in  w  at  400  mb.  should  not  seriously  affect  the  results. 


3030  700 mb. 

1420  850  mb. 

0  lOOOmb. 


1=12 


Fig.  27.  Vertical  cross-section  of  model  atmosphere, 
showing  lsobaric  surfaces  and  geopotential 
heights  used  in  testing  vortlcity  and  omega 
equation  methods  of  vertical  velocity 
calculations. 
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Fig.  31.  Same  as  Fig.  3°  except  for  400  mb. 


3.2,2  Vertical  Velocity,  Moisture  and  Cloud 

Vertical  velocities  were  calculated  for  19/0000  Z 
to  31/1200  Z,  July,  1973  at  the  inner  8x8  grldpoints 
of  Figure  3.  Cloud  observations,  classed  as  low, 
middle  of  clear,  were  then  taken  at  these  gridpoints 
from  the  satellite  images.  Also  noted  was  the  moist 
layer  thickness  observed  by  the  radiosonde  station 
nearest  to  a  given  gridpoint,  and  in  the  same  cloud 
or  clear  area. 

The  results  for  the  vorticity  method  of  u>  - 
calculations  are  shown  in  Tables  V-VIII.  For  Table  V, 
the  labels  low,  middle,  clear  1,  and  clear  2  are, 
respectively,  low  cloud  using  850-mb.  u;  -values, 
middle  cloud  using  700-mb.  UJ  -values,  and  clear 
using  850-mb.  and  700-mb.  uj  -values.  The  body  of 
the  table  contains  the  percentages  of  observations 
for  each  vertical  velocity-moist  layer  combination 
with  the  total  number  of  observations  shown  in 
brackets. 

Omega  equation  vertical  velocities,  calculated  at 
775  mb.  and  600  mb.,  are  shown  in  Tables  VIII  and  IX. 


Percentage  occurrences  of  combinations  of  vertical  velocity 
(by  vorticlty  method)  and  moist  layer  thickness  ranges 
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Table  VI 


Percentage  occurrences  of 
moist  layer  thickness  ranges 


Moist  Layer  Thickness 
( feet) 

Low 

(81) 

Middle 

(71) 

Clear 

(83) 

0£M*2000 

3.7 

8.5 

89.1 

20001M^8000 

67.9 

22.5 

10.8 

M*8000 

28.4 

69.0 

0.0 

Table  VII 


Percentage  occurrences  of  vertical  velocity  ranges 

using  vorticlty  method 


Vertical  Velocity 
(ub.sec.  ) 

Low 

(115) 

Middle 

(92)- 

Clear  1 
(97) 

Clear 

(97) 

V£-1.0 

0.0 

12.0 

1.0 

7.2 

-l.CKV£-0.5 

9.6 

15.2 

7.2 

16.4 

-0.5*V*0.0 

37.9 

26.1 

48.4 

36.0 

0.CKV^0.5 

46.1 

26.1 

39.1 

20.6 

0. 5^1.0 

7.0 

10.9 

2.0 

12.3 

V>1.0 

0.0 

9.8 

2.0 

7.2 

Same  as  Table  V  except  using  the  omega  equation  vertical  velocity  results 
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Table  IX 


Same  as  Table  VII  except 
using  omega  equation  vertical  velocity  results 


Vertical  Velocity 
(ub.sec.  ) 

Low 

(115) 

Middle 

(92). 

Clear  1 
(97) 

Clear 

(97) 

Vi  -1.0 

1.7 

5.4 

0.0 

0.0 

-l.o<v*-0.5 

5.2 

16.3 

1.0 

2.1 

-0.5^v*  0.0 

35.7 

45.7 

44.3 

41.2 

O.CUV^  0.5 

41.7 

22.8 

45.4 

45.4 

o.5iV*  l.o 

13.0 

9.8 

7.2 

7.2 

V>  1.0 

2.6 

0.0 

2.1 

4.1 

79 


For  the  vorttcity  method,  correlations  of  cloud 
with  moisture  and  ascent,  or  of  clear  areas  with  descent 
and  no  moisture  appear  poor.  For  example,  for  the  low 
cloud  category  In  Table  V  ,  In  the  2000  to  8000  feet 
range,  there  were  more  observations  of  weak  subsidence 
(35.8  percent)  than  weak  ascent  (19.8  percent).  For 
the  clear  2  cases,  in  the  0  to  2000  feet  range,  30.1 
percent  of  the  cases  show  weak  ascent,  and  only  20,4 
percent  show  weak  descent.  These  poor  correlations 
are  also  confirmed  In  Table  VII,  However,  in  Table  VI, 
it  may  be  seen  that  the  cloud  observations  are  discriminated 
much  better  by  moisture  alone.  That  is,  the  largest 
percentages  of  clear,  low  and  middle  cloud  cases  have 
moisture  thicknesses  0-2000  feet,  2000-8000  feet  and 
greater  than  8000  feet  respectively. 

From  the  omega  equation  results,  the  ascent-cloud 
and  descent-clear  correlations  are  improved.  In 
T&ble  VIII,  both  clear  categories  for  moisture  thickness 
0  to  2000  feet  show  more  subsidence  than  ascent.  For 
middle  cloud  with  M^  8000  feet,  more  cases  of  ascent 
are  observed.  However,  for  the  low  cloud  cases  for 
2000  £  M  ^  8000  feet,  more  subsidence  occurs.  Table  IX 
also  shows  this  improvement,  as  the  clear  areas  show 
more  subsidence  and  the  middle  cloud  areas  more  ascent, 
Ajrain,  the  low  cloud  category  has  few  cases  of  ascent. 
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Possible  reasons  for  this  slight  improvement  may  be 
the  followings  since  a  centred  time-derivative  of 
absolute  vorticity  Q  must  be  calculated  in  the  vorticity 
method  with  Q  calculated  every  twelve  hours,  an  error 
may  result  in  o>  due  to  24-hour  time-derivative  smoothing. 
Also,  the  thickness  term,  present  in  the  omega  equation, 
allows  for  thermal  advection  while  in  the  vorticity 
method,  no  thermal  advection  can  occur.  In  the  case 
study  in  Chapter  IV,  the  affect  of  this  term  is  shown. 

For  both  vertical  velocity  methods,  the  low  cloud 
results  are  the  worst.  This  may  be  expected  since 
775  mb.  and  850  mb. ,  the  levels  used  for  the  cloud-vertical 
velocity  calculations,  are  often  above  the  low  cloud 
and  thus  are  not  representative  levels  for  the  cloud. 


CHAPTER  TV 


CASE  STUDY  23/1200  Z,  JULY, 1973 

4.1  Introduction 

The  purpose  of  this  case  study  is  to  examine  cloud 
patterns  on  satellite  infrared  Images  and  correlate  them 
with  synoptic  features  such  as  lows,  troughs  and  TROWALS, 
as  well  as  with  the  distribution  of  moisture,  vorticity 
and  vertical  velocity. 

This  case  was  chosen  because  the  cloud  vortex, 
associated  with  a  500-mb.  low  centre,  was  located  over 
the  Arctic  Islands  where  conventional  data  coverage  was 
adequate  for  comparative  analysis  with  satellite  infrared 
Images  of  sufficient  quality  and  resolution.  The  case 
study  is  part  of  the  19/0000  Z  to  31/1200  Z,  July,  1973 
period,  for  which  geopotential  height  data  and  vertical 
velocities  were  available. 

4.2  Surface  Analysis 

The  surface  analysis,  Figure  32(b),  and  the  satellite 
images,  Figures  33,  34  show  a  cloud  vortex  southwest 
of  Thule  (THU) ,  close  to  the  surface  low  over  southern 
Ellesmere  Island.  The  TROWAL  over  northwestern 
Greenland  (see  Figure  32(b))  shows  up  as  a  band  of 
cloud  on  the  satellite  imaeres,  and  as  tongues  of  warm 
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air  on  the  850-  and  700-mb.  charts,  Figures  35(b)  and 
36(b)*  The  bright  white  cloud  over  the  Foxe  Basin  is 
associated  with  a  surface  wave  east  of  Hall  Beach  (YUX) 
and  a  1008-mb.  low  centred  near  Baker  Lake  (YBK). 

The  cold  front  extending  from  the  1008-mb.  low 
southwestward  is  well  marked  on  the  satellite  image 
by  a  line  of  clouds.  A  band  of  cloud  appears  also 
to  the  north  of  the  surface  warm  front,  over  the 
Beaufort  Sea. 

4.3  850-nib.  Analysis 

On  the  850-mb.  chart,  Figure  35(a),  areas  of  low 
cloud  are  outlined  by  wavy  lines  and  labelled  C(L). 
Moist  areas  with  dewpoint  depressions  DD^5*C  are 
labelled  M.  The  moist  areas  near  Mould  Bay  (YMD), 
northeast  from  Resolute  Bay  (YRB)  and  along  the 
Alaskan  coast  coincide  well  with  the  areas  of  low 
cloud.  The  low  cloud  from  Hall  Beach  to  north  of 
Great  Slave  Lake  and  northwestward  occurs  in  what 
appears  to  be  a  dry  area.  However,  since  this  is  an 
area  with  very  few  weather  observing  stations, 
correspondence  between  the  cloud  and  moisture  patterns 
is  likely  to  be  unsatisfactory. 

4.4  700-mb.  Analysis 

At  700  mb.,  Figure  36(a),  the  C(M )  areas  denote 

The  cloud  and  moisture  regions  coincide 


middle  cloud. 
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quite  well  over  Ellesmere  Island,  from  Baker  Lake 
southwestward,  and  in  the  Norman  Wells  ( YVQ)  area. 

However,  northeast  of  Norman  Wells  and  over  the  Beaufort- 
Sea,  middle  cloud  appears  on  the  satellite  image  but  no 
moisture  correlations  can  be  made  without  radiosonde 
stations  in  these  areas. 

4.5  500-mb.  Analysis 

A  strong  short-wave  trough,  from  the  vortlcity 

-5  -/ 

centre  22  x  10  sec.  southward  to  Hall  Beach,  has 
positive  vortlcity  advection  and  a  cloud  shield  to 
the  east  of  it.  A  second  trough,  from  west  of  Cambridge 
Bay  (YCB)  to  Great  Slave  Lake  has  a  weak  positive  vortlcity 
advection  area  and  narrow  bands  of  cloud  east  of  it. 
Upstream  from  this  trough,  there  is  negative  vortlcity 
advection  and  middle  cloud.  This  will  be  considered 
further  in  Section  4.7. 

4.6  Vertical  Velocity-Low  Level 

Since  low  cloud  frequently  is  below  the  775-rob. 
and  850-mb.  levels,  the  vertical  velocity  and  low  cloud 
patterns  may  not  correlate  well.  This  is  evident  in 
Figures  38(a)  and  39(a).  The  low  cloud-ascent  areas 
coincide  over  Northern  Greenland,  Ellesmere  Island  and 
from  Hall  Beach  to  north  of  Baker  Lake  but  not  in  the 
aroas  west  of  Baker  Lake  to  Great  Bear  Lake  and  Inuvlk. 


Fig.  37.  500-rob.  contour  and  vorticlty 

chart  for  23/1200  Z,  July,  1973. 
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Vertical  velocity  field  by  vorticity  method  Vertical  velocity  field  by  vorticity  me 
at  H50  mb.  and  low  cloud  for  23/1200  Z,  July,  at  700  mb.  and  middle  cloud  for  23/1200 
1973.  July,  1973. 
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Vertical  velocity  field  by  omega  equation  Vertical  velocity  field  by  omega  equation 
method  at  775  ®b.  and  low  cloud  for  method  at  600  mb,  and  middle  cloud  for 

23/1200  Z,  July,  1973.  23/1200  Z,  July,  1973. 


In  the  central  Arctic  Islands  the  omega  equation 
vertical  velocities  show  better  ascent-cloud  correlations. 

4,7  Vertical  Velocity  -  Middle  Level 

Figures  38(b)  and  39(b)  show  the  ?00-mb.  and  600-mb. 
vertical  velocity  patterns.  At  each  level,  middle  cloud 
and  ascent  occur  in  the  Ellesmere  Island-northern 
Greenland  areas  but  the  omega  equation  results  indicate 
closer  vertical  velocity-cloud  correlations.  Ascent 
and  middle  cloud  also  occur  together  from  Hall  Beach 
to  the  vicinity  of  Fort  Smith  (Y3M).  However,  both 
vertical  velocity  calculations  generate  subsidence  in 
the  vicinity  of  75°N  between  80°W  and  100°W  and  from 
Great  Slave  Lake  to  Inuvik.  The  major  difference 
between  the  vertical  velocity  calculations  appears 
north  of  the  Alaskan  coast.  At  500  mb.,  since  negative 
vorticity  advectlon  is  occurring,  subsidence  is  expected 
and  is  evident  in  Figure  38(b).  However,  the  omega 
equation  generates  ascent  in  this  same  area.  Since 
the  thermal  pattern  at  700  mb.  (see  Figure  36(b)) 
indicates  warm  air  advection  and  a  thickness  advectlon 
term  is  included  in  the  omega  equation,  it  is  this 
warm  air  advection  which  is  causing  the  ascent.  Thus, 
the  omega  equation  is  superior  in  this  instance  because 
it  includes  both  thickness  advection  and  vorticity. 


CHAPTER  V 


PROBLEMS  ENCOUNTERED  AND  SUGGESTIONS  FOR  IMPROVEMENTS 


5*1  Cloud  Classification 

As  noted  earlier  in  Chapter  I,  the  cloud  classification 
derived  from  satellite  infrared  images  is  based  on  a 
"grey-scale"  of  varying  shades  corresponding  to  the 
radiating  temperatures  of  the  cloud  masses  scanned. 
Occasionally  problems  may  arise  with  broken  middle  cloud, 
if  the  cloud  and  clear  patches  are  below  the  resolution 
of  the  Infrared  scanner.  According  to  Anderson  et  al  (1973) t 
if  there  are  several  sources  of  radiation  below  the 
resolution  of  the  sensor,  the  receiver  integrates  the 
radiances  into  a  single  response.  Thus,  for  broken 
middle  cloud,  the  scanner  will  respond  to  an  average 
temperature,  and  produce  an  Intermediate  shade  of  grey 
composed  of  the  dark  grey  of  the  warm  Earth* s  surface, 
and  the  white  of  the  cold  cloud  mass.  The  resulting  shade 
of  grey  may  lead  to  an  incorrect  classification  of  middle 
cloud  as  low  cloud.  Similar  problems  have  been  encountered 
by  Glahn  (1966).  He  found  that  individual  cloud  elements 
with  a  diameter  slightly  less  than  the  resolution  area  of 
the  infrared  scanner  can  have  their  tops  Judged  to  be 
several  thousand  feet  below  their  actual  height. 
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A  problem  also  arises  in  the  classification  of 
clouds  as  only  middle  or  low,  depending  on  their  crrey 
shade.  A  more  precise  determination  of  the  cloud  top 
heights  would  be  an  improvement  as  then  the  meteorological 
parameters  associated  with  the  cloud  could  be  found  more 
accurately.  A  possible  way  to  find  the  cloud  top  heights 
would  be  to  use  a  grey  scale  for  comparison  with  satellite 
image  grey  shades.  For  an  image, -the  grey  shade  of  a 
surface  at  a  known  temperature  (e.g.  an  open  water  area 
under  clear  skies)  could  be  compared  to  the  grey  scale. 

The  temperatures  corresponding  to  other  shades  of  grey 
on  the  scale  could  be  found  relative  to  this  known 
'‘reference"  grey  shade  temperature.  Then,  the  shades 
of  other  areas  on  the  satellite  image  could  be  compared 
to  the  grey  scale  to  determine  the  temperatures  of  the 
radiating  areas.  From  radiosonde  reports,  the  cloud 
top  heights  could  then  be  found  along  with  the  corresponding 
meteorological  parameters. 

5,2  Vertical  Velocity  -  Cloud  Correlations 

As  was  evident  in  the  vertical  velocity  contingency 
tables  given  in  Chapter  III,  the  expected  correlations 
between  ascent  and  low  cloud  were  often  poor.  One 
reason  for  this  is  that  the  vertical  velocity  is 

level  above  the  low  cloud  layer.  To 


calculated  at  a 
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improve  the  correlations,  vertical  velocities  produced 
by  surface  terrain  and  friction  In  the  boundary  layer 
should  be  calculated  and  correlated  with  the  low  cloud. 

Barr  et  al.  (1966)  studied  patterns  of  vertical 
velocity  associated  with  cloud  vortices  over  the 
northeastern  United  States,  They  found  that  both 
horizontal  and  vertical  motions  are  Important  in  large- 
scale  cloud  systems.  The  vertical  velocities,  calculated 
using  the  omega  equation,  show  reasonable  agreement  with 
cloud  patterns.  However,  the  authors  noted  that  since 
the  cloud  element  sizes  on  the  satellite  images  can  be 
several  orders  of  magnitude  smaller  than  the  scale  at 
which  the  vertical  velocities  are  calculated  (in  their 
study,  the  grid-length  was  165  km.),  only  the  broader 
aspects  of  the  cloud  patterns  may  be  expected  to  correlate 
well  with  the  vertical  motions.  In  this  study,  with  a 
grid-length  of  381  km. ,  the  point-by-point  comparisons  of 
vertical  velocities  and  cloud  showed  generally  poor  results. 
However,  in  the  case  study,  the  broad  aspects  of  the 
cloud  and  vertical  velocity  patterns  coincided  fairly 
well. 

Barr  et  al  (1966)  also  stressed  the  Importance  of 
horizontal  motions  in  the  maintenance  of  cloud  systems. 

In  their  analysis,  a  simple  trajectory  model  was  used  in 
which  air  oarcels  are  advected  by  the  wind  along  a  surface 
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of  constant  pressure.  In  this  way,  a  rough  estimate  of 
the  history  of  the  vertical  motion  of  the  parcel  may 
be  made  before  the  time  of  the  cloud  analysis.  This 
was  done  by  observing  the  vertical  motion  areas  for 
several  hours  before  the  time  of  the  cloud  observations. 
They  found  that,  frequently,  a  mass  of  cloud,  observed 
within  an  area  of  descent  at  the  time  of  the  vertical 
velocity  calculation,  may  have  been  in  an  area  of  ascent 
for  several  hours  prior  to  the  observation  time.  In 
such  cases,  the  net  vertical  velocity  can  be  ascent, 
while  the  descent,  which  may  be  occurring  at  the  time 
of  the  cloud  observation  is  too  short  or  too  weak  to 
dissipate  the  cloud. 

Barr  et  al  (1966)  also  found  that  the  vertical  motion 
correlated  best  with  the  cloud  patterns  early  in  the  vortex 
development  stage,  before  an  extensive  cyclonic  circulation 
developed.  Once  a  closed  circulation  has  been  set  up,  dry, 
clear  air  rotates  around  the  eastern  areas  of  the  vortex 
into  the  ascent  area,  and  moist,  cloud-bearing  air 
circulates  around  the  northern  parts  into  the  area  of 
descent  west  of  the  vortex  centre.  The  500-mb.  lows 
examined  in  the  present  study  had  closed  and  well-developed 
circulations  which  usually  persisted  for  several  days. 

Since  quasi— rotary  motion  of  cloudy  and  clear  areas 
within  a  vortex  is  a  common  feature  of  the  atmospheric 
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circulation,  the  use  of  even  a  fairly  crude  trajectory 
model  may  improve  the  correlations  between  the  vertical 
motion  and  the  cloud  patterns. 

A  final  point  concerns  the  nature  of  the  vertical 
motion  calculation.  The  superiority  of  the  omega 
equation  method  has  already  been  demonstrated  by  its 
ability  to  Include  thickness  advection.  However,  neither 
the  omega  nor  the  vorticity  method  include  diabatic 
terms  in  the  vertical  velocity  determinations.  Bullock 
et  al  (1969)  compared  vertical  motions  resulting  from 
diabatic  and  adiabatic  processes  separately.  The 
diabatic  process  considered  was  terrestrial  radiation. 

They  found  that,  for  clear  skies,  the  Infrared  radiation 
contributes  considerable  downward  motion  to  the  alrmass. 

Since  vertical  velocities,  determined  in  the  present 
study,  are  quite  small,  this  1  positive  vertical  velocity 
due  to  the  Infrared  components,  if  included,  could  have 
a  significant  effect  on  the  adiabatic  vertical  velocities. 

It  could  reverse  the  vertical  velocity  sign  from  negative 
to  positive  or  it  could  increase  an  already  positive 
vertical  velocity.  For  cloudy  skies,  the  mean  infrared 
component  of  the  vertical  velocity  in  the  Bullock  et  al  (1969) 
study  shows  negative  vertical  velocities  just  below  the 
cloud  base,  and  positive  values  just  above  the  cloud  layer. 
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Thus,  this  infrared  component  could  have  an  influence 
on  the  net  vertical  velocity,  depending  on  the  level 
at  which  it  is  calculated  relative  to  cloud  layers. 


CHAPTER  VI  . 


SUMMARY  AMD  CONCLUSIONS 


There  are  several  problems  in  Arctic  meteorology 
for  which  satellites  offer  assistance.  The  sparseness 
of  data  is  perhaps  the  foremost  difficulty.  This  is 
particularly  evident  over  the  Beaufort  Sea  and  the  vast 
Arctic  Ocean,  sources  of  many  low  pressure  areas  that 
track  across  the  Queen  Elizabeth  Islands.  Satellite 
coverage,  coupled  with  the  recognition  of  cloud  patterns 
and  associated  meteorological  features,  provide  for 
a  considerable  improvement  in  the  analysis  and  forecasting 
for  these  areas.  Satellite  image  analysis  is  also 
valuable  in  assessing  the  extent  of  cloud  masses  in 
more  data-rich  areas.  Fog  and  stratus  which  may  be 
reported  at  several  high  Arctic  stations,  could  be 
coastal,  widespread  or  topped  by  much  thicker  cloud. 
Satellite  information  helps  to  differentiate  among 
the  three  cases.  In  the  Arctic  winter  special  problems 
arise  where  satellite  photos  are  useful.  Very  low 
temperatures  make  radiosonde  soundings  of  moisture 
unreliable  for  the  determination  of  saturated  atmospheric 
layers.  Again  infrared  satellite  Images  are  helpful 
in  overcoming  this  problem  by  the  delineation  of 
wintertime  cloud  masses. 
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The  present  study  has  considered  four  aspects  of 
Arctic  meteorology: 

(1)  cyclone  tracks  at  500-mb. , 

(2)  moisture  and  temperature  distributions  in 
cloud  masses  associated  with  upper-level  cyclones, 

(3)  vertical  velocities  and  moisture  thicknesses 
in  cloud  areas  around  vortices,  and, 

(4)  the  correlation  of  synoptic  features  and 
cloud  areas  from  a  detailed  case  study. 

Figures  11  to  16  demonstrate  the  persistence  and 
complexity  of  the  500-mb.  cyclone  paths.  Occasionally, 
a  track  in  a  sparse  data  area  appears  complex,  but  the 
apparent  complexity  may  be  the  result  of  the  500-mb. 
objective  analysis  scheme.  Several  more  years  of  data 
would  be  required  to  determine  seasonal  mean  500-mb.  low 
positions  or  favoured  tracks  for  these  systems. 

To  determine  the  moisture  and  temperature  distributions 
associated  with  the  cloud  masses  spiralling  around  these 
cyclones,  a  crude  cloud-band  model  was  used  where  clear 
and  cloud  areas  around  a  vortex  centre  were  categorized. 
Subsequently,  radiosonde  reports  from  upper  air  stations 
located  in  the  various  areas  in  the  cloud  band  model 
were  used  to  obtain  moisture  and  temperature  distributions, 
as  well  as  central  tendencies  and  dispersion  statistics 
for  these  areas.  The  dewpoint  depression  distributions 
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were  as  expected,  that  Is  low  cloud  areas  showed  positively 
skew  distributions  at  the  surface  and  850-mb.  and  middle 
cloud  areas  had  positively  skew  distrubutions  at  700-mb. 

The  surface  and  850-mb.  level,  and  the  700-mb.  level 
were  used  as  rough  approximations  for  the  low  and  middle 
cloud  levels,  respectively.  Low  700-mb.  dewpoint 
depressions  were  occasionally  observed  in  low  cloud 
regions,  likely  caused  by  moisture  extending  above  these 
low  cloud  areas.  Also,  moisture  sometimes  extended  below 
the  middle  cloud  areas,  resulting  in  positively  skew 
surface  and  850-mb.  dewpoint  depression  distributions 
in  middle  cloud  regions. 

The  uncertainty  of  the  cloud  height  determined  from 
the  satellite  images  appears  to  be  a  major  problem 
encountered  in  this  study.  An  alternative,  using  a 
grey  scale  to  determine  the  cloud  top  temperatures  and 
hence  their  heights,  was  suggested  in  Chapter  V.  For 
future  studies,  this  should  be  a  considerable  improvement 
over  the  subjective  grey  shade  determination  of  low  or 

middle  cloud  as  used  here. 

To  calculate  vertical  velocities,  geopotential 
height  data  on  a  12  x  12  grid  were  used,  for  a  two  week 
period  in  July,  1973.  Contingency  tables  were  constructed 
to  determine  the  percentages  of  occurrence  of  combinations 
of  vertical  velocities  and  moist-layer  thicknesses.  i he 
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results  of  this  work  were  generally  not  promising  in 
that  regardless  of  the  moisture  layer  thicknesses,  a 
larger  percentage  of  subsidence  was  calculated  for  low 
cloud  areas  and  more  ascent  for  clear  areas,  the  opposite 
to  that  expected.  Middle  cloud  results  were  better,  with 
a  slightly  larger  percentage  of  cases  showing  ascent. 
Several  reasons  for  these  results  have  been  considered. 

The  use  of  the  vertical  velocities  at  850  and  775  mb.  , 
and  700  and  600  mb,  as  representative  of  low  and  middle 
cloud  levels,  respectively,  is  open  to  question.  Again, 
a  more  exact  estimate  of  the  level  of  the  cloud  must  be 
found  and  used  for  the  vertical  velocity  level. 

Computation  of  terrain- induced  vertical  velocities,  as 
representative  of  the  vertical  motion  in  low  cloud  areas, 
would  likely  increase  the  percentage  of  cases  of  ascent 
in  low  cloud  areas,  since  the  latter  are  frequently 
surface  based. 

The  methods  of  computation  of  the  vertical  motion 
have  also  been  questioned.  Both  calculations  include 
only  adiabatic  processes  of  vertical  motion,  ignoring 
diabatic  processes  such  as  infrared  radiation,  which, 
in  clear  sky  conditions,  contributes  a  positive  vertical 
velocity  component.  In  cloud  areas  the  infrared  component 
has  a  mean  negative  value  just  below  the  cloud  base  and 
a  mean  positive  value  just  above  it..  Thus,  depending  on 
the  level  of  computation  of  the  vertical  velocity,  the 
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Infrared  component  could  have  an  effect  in  cloudy  conditions. 

Finally,  an  air  parcel  trajectory  determination  has 
been  suggested.  Since  500-mb.  cyclones  have  frequently 
spiralling  bands  of  cloud  associated  with  them,  dry, 
clear  air  may  be  intruding  into  ascent  regions,  and  moist, 
cloud  air  into  subsidence  regions  at  the  time  of 
observations  of  the  satellite  image.  A  trajectory 
model  could  be  used  to  determine  an  approximation  of 
the  net  vertical  motion  an  air  parcel  has  experienced 
prior  to  the  satellite  image  time. 

A  detailed  examination  of  a  case  has  been  done  to 
correlate  cloud  patterns  with  synoptic  features.  In 
regions  with  upper-air  data,  the  700-mb.  and  850-mb. 
moist  areas  coincided  fairly  well  with  the  satellite 
cloud  areas.  At  500  mb. ,  the  positive  vorticlty  advection 
areas,  in  most  cases,  had  cloud  associated  but,  in  one 
instance,  negative  vorticlty  advection  and  cloud  occurred 
together.  For  this  same  area,  positive  vertical  velocities 
were  calculated  by  the  vorticlty  method,  whereas  the 
omega  equation  calculations  correctly  indicated  negative 
vertical  velocities.  This  superiority  of  the  omega 
equation  is  partly  the  result  of  the  thickness  advection 
term.  Since  warm-air  advection  was  also  occurring  in 
the  negative  vorticlty  advection  area  in  question,  the 
omega  equation  correctly  predicted  ascending  motion. 


. 
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It  Is  hoped  that  the  results  of  this  study  will  be 
useful  for  guidance  as  a  first  approximation  of  the 
meteorological  parameters  associated  with  500-nib.  cyclones. 
Also,  perhaps  the  problems  and  pitfalls  discussed  will 
be  of  assistance  to  others  undertaking  similar  studies 
using  satellite  Images  in  the  high  Arctic. 
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APPENDIX  A 


A  "LOG -NORMAL.!  TRANSFORMATION 


The  transformation  is  for  the  surface  dewpoint 
depression  distribution  for  the  C(L)  area  shown  in 
Figure  17. 

(a)  The  skewness  ft  and  kurtosis  ft  are  calculated 
as  follows: 

For  the  C(L)  surface  dewpoint  depression  distribution, 
the  statistics  are  given  below. 

Mean  (x)  =  2.9 
Variance  {juz  )  =  7.2 
Standard  Deviation  (05)  =  2.7 
Third  moment  {  M3)  -  58.0 
Fourth  moment  (^u¥)  =  710.6 
From  equation  (1.27)* 
ft  *  +  3  •  0 

From  equation  (1.28) 

ft*+10.6 

Since  both  ft  and  ft*  are  positive,  the  "log-normal” 
transformation  may  be  applied. 

(b)  Using  the  rx  given  above,  tu  is  calculated 
from  equation  (1.3D* 

w=  0. 8 

(c)  <TU  is  calculated  from  equation  (1.32). 

<r*=  0.7 
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(d)  Calculate from  equation  (1.30)  for  each 
value  of  the  Independent  variable  x(here  the  dewpoint 
depression) . 

From  a  table  of  probability  integrals  for  the 
normal  frequency  curve  such  as  Appendix  11(c)  in 
Brooks  et  al  (1953) »  the  frequency  distribution  of 

was  obtained.  The  original  dewpoint  depression 
frequencies  and  the  transformed  variable  frequencies 

are  shown  in  Table  A-I  and  plotted  in  Figure  26. 
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Table  A.- 1 

Dewpoint  depression  and 
transformed  variable  frequencies 


Dewpoint 

Frequency 

Transformed  Variable 

Frequency  of 

Depression 

Transformed 

(°c.) 

Variable 

0.0 

6 

-2.6  • 

18.6 

1.0 

25 

-0.9 

26.0 

2.0 

23 

-0.1 

18.8 

3.0 

20 

0.4 

11.9 

4.0 

5 

0.8 

7.6 

5.0 

9 

1.0 

4.8 

6.0 

3 

1.3 

3.0 

7.0 

1 

1.5 

2.1 

8.0 

3 

1.7 

1.4 

9.0 

0 

1.8 

0.9 

10.0 

0 

2.0 

0.7 

11.0 

0 

2.1 

0.6 

12.0 

0 

2.2 

0.4 

13.0 

1 

2.3 

0.3 

14.0 

1 

2.4 

0.2 

15.0 

1 

2.5 

0.2 

